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PART I 


CONSTRUCTION OF THE CHRONIC 2G CENTRIFUGE 


In formulating design characteristics of the Chronic 2g Centrifuge, 
several constraints were considered in order to ensure the constancy and 
reliability of vestibular data retrieved from the rats housed on board 
the centrifuge. A prime concern was that of reliability. The data 
would have little meaning if the unit were inoperable for any amount of 
time, requiring repair or maintenance. It was predetermined that any 
maintenance required should be scheduled maintenance which could be 
administered at intervals of at least several months. The centrifuge 
must also run continuously and at a constant velocity. The centrifuge 
was to have as large a radius arm as possible in order to minimize the 
coriolis phenomena effects that could alter vestibular system response. 
Both the velocity and radius, however, were to be infinitely variable, 
within certain specified limits, in order to provide flexibility in 
experimental design if the need arose at a future date. The centrifuge 
was to provide living space for a minimum of 20 rats, which were to be 
housed in standard 18" x 10" x 7" deep laboratory cages. The rats were 
to receive food and water ad libitum . It was emphasized that wobble 
and vibration be kept to a minimum. The cages were to be oriented on 
the radius arm such that at any speed, from zero to maximum velocity, 
the floor of the cage was to be perpendicular to the resultant gravity 
vector . 

The Electrical Shop in the Department of Electrical Engineering at 
The Ohio State University agreed to design and construct the centrifuge. 
They agreed to a six-month deadline but requested that the design stage 
be kept to a minimum in order to assure completion with the six-month 
period. This request proved to be an asset in speeding construction and 
completion of the centrifuge within the allotted time span. 

The dimensions of the centrifuge were dictated by the size of the 
roam in which it was to be housed. It was decided that the maximum 
radius was to be 10.5 ft. By dividing the length of the circumference 
by the length of each cage, it was determined that we could easily 
accommodate 20 cages at the maximum radius . At shorter radii fewer 
cages can be accommodated at any one particular radius because of the 
reduced circumference. We have the option, however, of staggering cages 
at several different radii. This will be available in later exper imen ts 
in which we wish to expose different populations of experimental animals 
to different g levels while maintaining a constant velocity of rotation 
and coriolis phenomena. The design includes the potential of positioning 
a single cage or a vertical stack of several cages over the center axis 
of the centrifuge, thus positioning the rats at essentially zero radius, 
lg. However, since they will turn at the same velocity as the rest of 
the centrifuge they will be exposed to the identical potential coriolis 
phenomena . 
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The standard polyethylene laboratory cages are individually held 
in an aluminum cage holder. The cage holder is formed and heliarc 
welded on a jig (Fig. l). The individual cage slips into the cage 
holder (Fig. 2) and can be easily removed for cage cleaning. This size 
cage can accommodate up to four adult rats or one rat with a litter. 

A gimbal cage support suspension device was accepted following 
tests on a prototype (Fig. 2). The gimbal consists of a ball and socket 
joint and a horizontal axis (Fig. 3). The ball is made of brass and 
can rotate sufficiently within the socket assembly to maintain alignment 
of the horizontal axis . The cage holder is suspended from the ball and 
socket by a short support stud riveted by means of a flange to the cage 
support bracket (Fig. 3 )- The stud slips through a hole in the brass 
ball and the assembly is secured by a circlip at the end of the stud. 

This gimbal has proved to present a minimum of friction and is satis- 
factory in maintaining the floor of the cage perpendicular to the 
gravity vector, regardless of speed of rotation (Fig. 4 ). 

Food and water are provided for the rats by a conventional pellet 
tray and water bottle holder incorporated into the cage lid. Inadver- 
tent water loss from the water bottle has not proved to be a problem. 

This may be partially due to the fact that we have selected dropper 
spouts with a relatively small diameter opening and have cut the spout 
shaft to keep it as short as possible and still allow the rat easy 
access to the water. 

Each cage support is anchored to its pair of radius arms or radius 
arm extensions by means of two pinch clamps, one directly below each 
gimbal. The clamp is secured by tightening a single bolt. The clamp 
holds the cage support in its position by friction and the system is 
designed so that the cage may be positioned at any radius from the hub 
to the outer limit of the radius am extensions . The clamps are aligned 
on the box channels making up the radius arms and radius arm extensions 
so that if both cage clamps should fail, which is unlikely, the cage 
would slide to the maximum radius of the radius arms or radius extension 
arms . They would not lose their attachment to the centrifuge and thus 
could not fly off the centrifuge producing a potential hazard to personnel. 

The simplest and most rigid radius arm construction was employed. 

Box aluminum alloy channel (2" x 4 ") was chosen for radius arms and 
radius arm extensions. This construction gives the greatest rigidity 
and strength per pound weight and, at the same time, simplifies the 
method of assembly. At critical points such as the point at which each 
radius arm is joined to the hub and where the extension arm is attached, 
a sleeve fit was used for strength and ease of assembly. The internal 
and external sections of each sleeve were then secured with bolts (Fig. 5). 

Each radius arm assembly was divided into two parts. The inner is 
the radius arm proper and the outer is the radius arm extension. This 
method was devised for ease of construction and also to provide more 


2 




Fig. 1 - Jig developed to form and weld aluminum cage 
holder 



Fig. 2 - Cage holder for standard laboratory- 
cage measuring l8" x 10" x 7" deep 
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Fig. 3 - Detail of gimbal unit, with brass ball and 
support stud, which allows swivel of cage 
holder on centrifuge 



Fig. 4 - Cages on test device rotating at about 100 rpm 
(Note cage tilt in response to combined gravi- 
tational and centrifugal forces) 
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flexibility in the operation of the centrifuge. With just the radius 
arms secured to the hub, the maximum effective radius of the centrifuge 
xs 5.25 ft (Fig. 6 ). This considerably reduces the weight of the 
centrifuge wheel and also reduces wind-loading. The power requirements 
for driving this effectively smaller centrifuge are reduced and thus, 
with the power unit available, higher rotational velocities can be 
achieved. The increase in rotational velocity could produce higher g 
levels or higher potential coriolis stimulation or both, on the experi- 
mental design. 

The radius arm extension is a "Y" shaped extension arm. The base 
of the "Y" plugs into the terminal end of the radius arm and extends 
the radius of the centrifuge to 10. 5 ft (Fig. 7). The use of the radius 
arm extension increases the total weight and wind- loading force of the 
centrifuge -wheel. The increased radius, however, reduces the potential 
coriolis effect generated at any particular g level produced. This has 
proved to be an important stimulus parameter when the vestibular system 
is being examined. 

The hub assembly is the critical part of the centrifuge. It is 
required to anchor the radius arms, transmit torque from the vertical 
shaft to the radius arms, and transmit the entire weight of the super- 
structure to the vertical drive shaft and its support bearing. In order 
to ensure success of the hub assembly construction a careful design was 
completed and a wooden mock-up was built (Figs. 8 and 9 ). It was 
decided that the hub be constructed as a single welded unit in order to 
provide strength. The same aluminum alloy was used in order to help 
reduce weight. The hub assembly is constructed of two 0.5-inch-thick 
plates 22 inches in diameter. They fit around a 4-inch outside diameter 
center sleeve 13.5 inches long. A steel quick-dis connect pulley hub 
fits into each end of the sleeve. These pulley hubs are keyed for a 
two-inch shaft. The circular plates are strengthened by five triangular 
gussets welded between the plate and center sleeve on both the upper 
and lower side of the assembly (Fig. 10). 

Twenty specially welded and machined box studs were incorporated 
into the hub assembly between the circular disks . Each stud is 15 inches 
long and is individually fitted to accept its corresponding radius arm, 
with a minimum of clearance (Fig. 10). The individual components of 
the hub are heliare welded to form an integral assembly. 

The base is made of four- inch steel channel welded into a six- 
foot square box frame (Fig. 11). Two parallel cross braces directly 
support the load and transfer it to the box frame dir ectly and also by 
means of four cross angle braces. In order to achieve added support 
and stability, each side of the box frame is provided with a bolt-on 
leg that adds ten feet to the support base in each direction (Fig. 12). 

The pedestal assembly incorporates the drive motor, speed reducer, 
drive shaft, and support bearings (Fig. 13). Four hollow steel cylinders. 
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Fig. 8 - Wooden mock-up of centrifuge hub with top 
support disk and gussets in place 
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Fig. 9 - Wooden mock-up of centrifuge hub with top 
support disk removed revealing position of 
studs 






Fig. 11 - Six-foot square centrifuge support base 

(Center post represents axis of centrifuge) 



Fig. 12 - Six-foot square centrifuge support base with 
cross braces and bolt-on legs 
(At left is pedestal assembly) 
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Fig. 14 - Completed centrifuge with cages holders and 
cages positioned at maximum radius 


10 




measuring 3-5 x 22.5 inches, support the base plates of the pedestal 
assembly. Four hollow steel cylinders, 3-5 x 9*5 inches, are bolted to 
the base plate and support the upper plate of the pedestal assembly. 

The steel base plate is 20 x 20 x 1 inch thick and is machined on the 
upper and lower surfaces to provide a smooth flat surface. The steel 
upper plate is 16 x 16 x 1 inch thick. 

The centrifuge wheel weight and maximum diameter were taken into 
consideration in estimating the power required to drive the centrifuge. 

It was found that the leading edge of all cages as they rotated displaced 
a considerable quantity of air, thus developing a wind-loading force 
opposing the rotation of the centrifuge. Several tests were attempted 
to determine a value for this force, and a rough estimate was determined. 
In addition, it was found that the wind-loading effect on the leading 
edge of the combined radius arms also contributed to this force. Based 
on these estimates it was concluded that a three-horsepower direct 
current motor would be adequate to drive the centrifuge within the 
speed ranges necessary to generate at least 2g at any point along the 
radius arm of the centrifuge. 

The drive motor selected is a General Electric, Type L186ATC, 
three-horsepower direct current motor. It is suspended in a vertical 
position in the base of the pedestal. It is coupled directly to a 
Boston Gear Optimount 17:1 flanged speed reducer. The speed reducer is 
bolted to the underside of the base plate and accepts the 2 7/l6-inch 
keyed drive shaft. The drive shaft reduces to a two-inch diameter above 
the bearing supports. The maximum motor speed is 1750 rpm. The speed 
reducer allows the output drive shaft to turn at a maximum of about 
100 rpm. Tests of the completed centrifuge have shown that maximum 
speed can be approximated with the radius arm and inner ring assembled 
(Fig. 14) . With the radius arm extensions added, the maximum speed is 
approximately 30 rpm. The only maintenance required is to change the 
speed reducer oil once every 40 days, and to periodically inspect and 
replace the motor brushes . 

Centrifuge speed is controlled by a Boston Gear Ratiotrol, produced 
by the Boston Gear Division of North American Rockwell. This device is 
a full-wave, solid-state SCR control and utilizes modular construction 
and transistorized circuitry, including integrated circuits. The control 
will accelerate the centrifuge from zero to a present velocity. It will 
maintain this speed continuously within ±1% of motor base speed. 
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PART II 


INTRODUCTION 


The vestibular system of equilibrium includes both afferent and 
efferent components. The sensory end organ is made up of two functional 
subunits. One of these is sensitive to gravity and linear acceleration 
and the second of these is sensitive to turning movements of the head, 
viz., angular acceleration. Linear acceleration is detected by the macula 
of the utriculus and possibly the macula of the sacculus. Angular accer- 
eration is detected by the cristae ampullares of the three semicircular 
canals. The vestibular end organ provides the central nervous system with 
information as to the position of the head with respect to gravity and the 
acceleration of the head in space. The efferent components of this system 
are related to antigravity musculature, muscles associated with movements 
of the head and neck, and extraocular muscles. In addition, the vestibular 
system has been associated with the autonomic nervous system and arousal. 

Two other sensory systems also play a role in maintenance of equili- 
brium. First, the visual system functions primarily to maintain a visual 
horizon. Second, the proprioceptive system provides information as to 
the position of the head with respect to the body, the position of limbs 
with respect to the body, and the position of the limbs with respect to 
the substratusm. All three of these systems function together in order 
to assess the position of the head in space, the position of the head 
with respect to the body, and the position of the limbs relative to the 
body and the substratusm. 

The exact role played by each of these systems, or by the individual 
parts of each system, is not yet clearly known. In this study an attempt 
was made to establish a quantitative measurement of the ability of the 
white rat to utilize these equilibrium systems in order to maintain his 
balance during a dynamic situation. 


MATERIALS AND METHODS 

Two groups of rats were utilized. One group, the lg controls, were 
adult rats that had been raised under standard laboratory conditions under 
the effects of the earth standard gravity of lg. The second group, the 
2g experimentals , were adult rats that had been maintained on a centrifuge 
for a minimum of 60 days, 24 hours per day. 

The centrifuge was a small table top device measuring about 1 m in 
diameter. In order to produce the required 2g environment, the centri- 
fuge was rotated at 60 rpm. The cages were tilted inward about 30° in 
order to position the 2g gravity vector perpendicular to the floor of the 
cage. Both groups of rats received food and water ad libidum. 
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The device used to test the dynamic equilibrium characteristics of 
the rats was designed after the rail device for monkeys, described by 
Igarashi (1968). The rat rail device measures 3 ft. long, If- ft. wide 
and 2 ft. high. Extending beyond the test box are two start boxes, one 
located at each end. Each is separated from the test box by a sliding 
door. A rail, 1 inch in diameter, runs the length of the box and is 
cross-shaped in cross-secrion. The floor of the start boxes and the 
test box is an electric shock grid. The rail is at the same level as 
the floor of the start boxes, and is 4 inches above the electrified grid 
floor of the test box. The rail is connected to a D.C. motor and the 
rate of rotation of the rail can be controlled from 0-25 rpm. 


PROCEDURE 

During a single testing session, a rat was placed in the start box 
for 15 seconds at one end of the rail. The sliding door was raised and 
a light turned on at the opposite end of the test box. After five 
seconds the grid in the start box was electrified. The rat was allowed 
to remain in the test box for up to 60 seconds. There were ten lg con- 
trol and ten 2 g experimental rats. Each rat was tested first in the 
stationary rail test; i.e., each animal was run for ten trials each day 
for ten consecutive days. Following the stationary rail test, each rat 
was then tested in the rotating rail test. Each animal was run for one 
trial at each of the following rail rotation speeds: 0, 5> 10, 15» 20, 

and 25 rpm. The procedure for the rotating rail test was continued for 
three consecutive days . All 2 g rats were returned to the 2 g centrifuge 
following completion of tests each day. All testing was done at lg. 

Three sets of data were collected from each group of animals for 
each trial. First, it was determined if the rat actively avoided the 
shock in the start box that was initiated five seconds after the sliding 
door was raised. Second, the total time spent on the rail was recorded. 
Third, it was determined if the animal remained on the rail during the 
entire one-minute trial. The last two sets of data are similar in that 
they both are a measure of the animal's ability to maintain balance. 


RESULTS 


Stationary Rail Test 


There was no significant difference between the lg and 2 g rats in 
their ability to balance using time spent on rail as a measure. Table 
I indicates that by the fourth day of testing, both groups spent the 
entire sixty seconds on the rail. During the first three days of testing 
the 2g animals tended to fall earlier in the trial, but the group dif- 
ferences are not significant. 
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When the groups were compared as to whether each rat spent the 
entire sixty seconds on the rail, significant differences were seen on 
the first two days of testing (day 1, X 2 = 82, p < .01; day 2, X 2 = 12.5, 
p < .01). On the first two days rats in both grqpps performed better at 
maintaining balance on the stationary bar. During these first two days, 
the 2g group appeared to improve at a greater rate. On the third and 
subsequent days, there was no significant difference between the two 
groups . 

When avoidance was examined it was found that the two groups differed 
significantly on all days but the first and fifth days. On all days, 
including the first and fifth days, the 2g rats avoided the shock more 
often than the 1 g rats. Both groups improved avoidance scores from 
the performance on the first day. The lg group appeared to reach a 
plateau on the fourth day while the 2g group appeared to reach a plateau 
on the second day. 


Rotating Rail Test 


The lg and 2g rats differed significantly in their performance on 
the rotating rail. When the time spent on the rail was compared for each 
group, it was seen that on each of the three days of testing, the lg 
group remained on the rail significantly longer (Table IV). When the 
performance was broken down into categories based on the rate of rotation 
of the rail., it was found that at faster rates of rail rotation, the 
differences between the two groups becomes more apparent (Table V). At 
rail rotation rates of 0 and 5 rpm, summed over the three days of test- 
ing, there is no significant difference between the lg and 2g rats. At 
10 rpm, summed over three days, the two groups are significantly differ- 
ent at the 95% level of confidence, while at 15, 20, and 25 rpm rates the 
groups are significantly different at the 99% level of confidence. 

When avoidance was examined, it was found that the two groups 
differed significantly on all three days of testing (Table VI ). The 
2g group consistently avoided more frequently than the lg group. 


DISCUSSION 

Comparing the differences seen in the lg and 2g groups, it appears 
that the rotating rail test is a more sensitive tool in differentiating 
between the two groups. No significant difference was seen between the 
two groups in the Stationary Rail Test even though there was a consistent 
numerical difference when time on rail 'was used as the discriminative 
parameter (Table I). Table II, which displays the number of trials when 
each rat spent the entire time on the rail, revealed a significant dif- 
ference between the two groups on the first two days . This is reflected 
in Table V at the 0 and 5 rpm rail rotation speeds, but when the rail 
speed increased above 5 rpm, the differences between the two groups 
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becomes more apparent. At the higher rail rotation speeds the differences 
between the two groups appears to be progressively greater as suggested 
by the increasing "+" values (Table V). This difference may be more 
significant at rail rotation speeds above 25 rpm.* 

The lg and 2g rats differed in their ability to balance. This was 
consistently obvious to the experimenter, even in the stationary rail 
test. The 2g rats appeared to have much difficulty in maintaining 
themselves on the rail. They appeared to "lose their balance" often 
but soon learned to stay on the rail even though they appeared to have 
much more difficulty than the lg animals. These subjective evaluations 
on my part were not reflected in the scores seen in the Stationary Rail 
Test. At higher rail rotation speeds in the Rotating Rail Test, however, 
these differences began to become apparent in the test scores . 

The data measuring time on rail, for the first three days in the 
Stationary Rail Test revealed a consistent, but not significantly dif- 
ferent, difference between the lg and 2g rats (Table I). Similar 
differences were seen for the first five days when the two groups were 
compared as to the average number of trials in which each rat spent the 
entire time on the rail. The data were significantly different for the 
first two days (Table II). The Rotating Rail Test revealed that the 
two groups did not differ at 0 and 5 rpm. But when the test became 
more demanding, i.e., at higher rail rotation speeds, the differences 
between the two groups became significant. 

The data from both tests reveal that both the lg rats and 2g rats 
improved their ability to stay on the rail with experience. This is 
especially apparent when the scores are compared in the Stationary Rail 
Test which ran for a longer time period. If the "+" values are compared 
in Table IV, it appears that both groups are improving performance at 
about the same rate. There is no reason, based on these data, to suggest 
that the 2g rats will approach the lg rats in their balancing ability. 

This raises two questions. First, given enough training (e.g., several 
weeks) will the 2g rats approach the performance level of the lg rats? 
Second, if the 2g rats were to remain in the lg environment, would they 
eventually approach the performance of the lg rats? How long would this 
take? 


It is premature at this time to speculate on the anatomical and/or 
physiological basis for the behavioral differences seen. The data might 
suggest that it is a system that is called upon to function more in a 
dynamic situation (viz,, semicircular canals). This maybe the case. 
However, even when balancing on the stationary rail, the 2g rat experi- 
ences more difficulty in maintaining balance, implying that other sen- 
sory mechanisms may be involved. 

The differences seen in avoidance behavior were unexpected. The 
2g group avoided consistently more frequently than did the lg group. 

Both groups improved their avoidance scores with experience (Table III). 
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The X 2 values remained high even at the end of ten days of testing, 
suggesting that the difference between the two groups was remaining 
constant. The cause for these differences is unknown. We might 
speculate that the 2g environment has imposed a stress factor that has 
altered the arousal level in the 2g rats. Possibly the answer may lie 
with the reported vestibular fibers projecting to the brain stem re- 
ticular formation, and the known influence of vestibular stimulation on 
visceral efferents. 



Table 1 - Stationary Rail Test, (lg rats and 2 g rats compared as to 

average time, in seconds, spent on rail during one 60 second 
trial. Ten rats per group were exposed to ten trials each 
day for ten days. The groups are not significantly different 
for any day) . 


Day 


Group 

1 

2 

3 

4 

5 6 

7 

8 9 

10 

lg 

49.0 

52.2 

60.0 

60.0 

60.0 60.0 

60.0 

60.0 60.0 

60.0 

2g 

38.0 

47.3 

56.3 

60.0 

60.0 60.0 

60.0 

60.0 60.0 

60.0 

t value 

1.46 

0.62 

4.68 






Table 2 _ 

Stationary Rail Test, (lg rats and 2 g rats compared as to 
average number of trials in which each rat spent entire 
time on rail. Ten rats per group were exposed to ten 
trials each day for ten days). 






Day 




Group 

1 

2 

3 

4 

5 6 

7 

8 9 

10 

lg 

8.2 

9.0 

9-6 

10.0 

10.0 10.0 

10.0 

10.0 10.0 

10.0 

2g 

1.8 

7-0 

9.8 

9.8 

9.8 10.0 

10.0 

10.0 10.0 

10.0 

X s 

82.0** 

12.5** 

0.57 







Table 3 - Stationary Rail Test, (lg rats and 2 g rats compared as 
to avoidance responses made each day. Ten rats per 
group were exposed to ten trials per day for ten days). 


Group 





Day 





1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

lg 

.2 

.2 

.6 

3.4 

1.8 

2.2 

2.6 

1.6 

2.4 

2.8 

2g 

.6 

3-6 

4.8 

5.2 

2.8 

4.6 

5-4 

4.6 

4.4 

5.0 

X s 

2.1** 

37 - 4 ** 

44 . 8 ** 

6.6* 

2.8 

12.8** 

16.1** 

21.2** 

8.8** 

10.2** 
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Table 4 - Rotating Rail Test, (lg rats and 2g rats compared as to 
average time, in seconds, spent on rail daring one 60 
second trial. Ten rats per group were exposed to ten 
trials each day for three days . Data are averaged for 
rail speed from zero to 25 rpm) . 


Group 


Day 


1 

2 

3 

lg 

45.2 

48.6 

52.7 

2g 

30.5 

36.5 

35-6 

t value 

3.49* 

3.75* 

3.38* 


Table 5 - Rotating Rail Test, (lg rats and 2g rats compared as to 
average time, in seconds, spent on rail during one 60 
second trial. Ten rats per group were exposed to ten 
trials each day for three days . Data are averaged for 
three days ) . 


Group 



Rail Speed (rpm) 



0 

5 

10 

15 

20 

25 

lg 

58.0 

54.1 

50.9 

48.9 

40.9 

39.2 

2g 

57-1 

47.8 

35.1 

29.9 

22.3 

19.8 

t value 

0.31 

1-33 

2.48*-* 

3.52* 

3.03* 

3.62* 


Table 6 - Rotating Rail Test, (lg and 2g rats compared as to 
average number of avoidance responses made each day. 
Ten rats per group were exposed to ten trials per 
day for three days). 


Group 


Day 


1 

2 

3 

lg 

2.2 

2.0 

2.8 

2g 

5.0 

3.5 

3.2 

X 2 

27.24* 

7.56* 

6.54* 


* significantly different at the 0.05 level | This reference applies to 
significantly different at the 0.01 level j Tables II through VI. 
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PART III 


INTRODUCTION 

PHYSIOLOGICAL REVIEW AND MATHEMATICAL BASIS 


1. Vestibular Apparatus 

The membranous labyrinth is located in the inner ear. It con- 
sists of the cochlear duct, three semicircular canals, and two relatively 
large sacs known as the utricle and the saccule. The cochlear duct is 
concerned with hearing and has nothing to do with equilibrium; the role 
of saccule in orientation and postural control is not clearly understood 
in higher animals. However, the utricle and the semicircular canals are 
important for maintaining equilibrium. 

a. Utricle — The macula is the sensory epithelium of the 
utricle and consists of sensory cells, a gelatinous substance and many 
small calcium carbonate crystals called otoliths . The otoliths are more 
dense than the surrounding structures and fluid and, thus, are differ- 
entially attracted by gravity. The otoliths are embedded in the upper 
surface of the gelatinous substance; below the gelatinous substance are 
located many sensory hair cells . Two types of sensory cells have been 
identified and are referred to as Type I and Type II. Small sensory 
cilia project from the sensory cells into the gelatinous substance. 

It is presumed that a shifting of the otoliths results in a shearing 
action imparted to the cilia of the sensory cells which action results 
in depolarization of the sensory cell. 

b. Semicircular Canals — The three semicircular canals, known 
as the anterior, posterior, and horizontal (or lateral) semicircular 
canals, are arranged at approximately right angles to each other, so 
that they roughly form an orthogonal system. The lateral canals lie 

in a plane tilted approximately 30 degrees above the horizontal. The 
two vertical canals lie approximately vertical to the plane of the 
horizontal canal and 90 degrees to each other. 

Each canal is in communication with the cavity of the utricle. 

At one end of each canal the duct dilates into an enlarged cavity, the 
ampulla. The neurosensory epithelium of each semicircular canal lies 
within the ampulla. 

The neurosensory structure in the ampulla is called the crista 
ampullaris . The crista is composed of an overlying gelatinous cupula 
and sensory hair cells below. The cupula occludes the ampulla. Cilia 
project from the sensory hair cells into the cupula. The cilia in each 
crista ampullaris are uniformly directionally oriented. 
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When the head is rotated in any direction, the endolymph in the 
semicircular canals tends to remain stationary due to inertia. The 
inertia force causes deflection of the cupula. Deflection of the 
cupula results in a shearing force transmitted to the cilia of the 
sensory cells. This action results in depolarization or hyperpolari- 
zation of the cell, depending upon the direction of cupular deflection. 
The appropriate neural signals are transmitted via the eighth cranial 
nerve to the brainstem, cause a sensation of rotation and induce the 
compensatory eye movement called nystagmus . 

c. Nystagmus --Nystagmus is the reflex eye movement due to a 
deflection of the cupula. When the head is accelerated about an axis 
projected vertically through the center of the head, the eyes move 
slowly in a direction opposite to the direction of rotation with a 
velocity equal to that with which the head is turning. This is the 
slow phase of nystagmus . This movement tends to fixate a visual image 
on the retina, and minimize the displacement of the image on the retina. 
The rapid return of the eye to a new starting point in the orbit is 
called the fast phase of nystagmus. 


2. System Analysis 

In the field of control systems analysis, one of the most 
important things is to identify the transfer-function of the system, 
which gives information about the functional organization of the mecha- 
nism. From the point of view of system analysis, we are dealing with 
a very complicated system which has probably a number of feedback loops 
in the central nervous system and cross-couples with other systems 
(i.e., otoliths) through the C.N.S. But, on a biological basis it is 
not possible to determine the open-loop transfer functions in great 
detail. We can develop an over-all model of the system which explains 
behavioral "input -out put" data and is consistent with known physiolog- 
ical structures . 

a. System Transfer Function — The simplified horizontal vesti- 

bulo-ocular reflex arc is shown in Fig. l(a). For convenience in 
analysis, the system can be grouped and arranged in several block dia- 
grams which is summarized in Fig. 1 (b). This system has multivariable 
inputs and a single output, but for easy manipulation we have to re- 
strict attention to one input and one output which is shown in Fig. 1 
(c). The relation between input and output is expressed by the transfer 
H(s) in the complex frequency domain: 



Gravitational 



Fig. 1 (a) The simplified horizontal vestibulo- ocular reflex arc in block diagram. 

(b) The arrows indicate the direction of information flow. 

(c) The simplified transfer function with input [cn(s)] and output [V(s)]. 












where 


$(s) = Laplace transform of the input angular acceleration, and 
V(s) = Laplace transform of the velocity of slew-phase nystagmus. 

Since the time function of the velocity of slow-phase nystagmus 
following an impulse input has the form of an exponential decay, the 
system transfer function H(s) can be assumed to be: 


H(s) 


b n _ 1 s n ~ 1 + b n _ 2 s n “ 2 + ... + b 2 s 2 + b s + b 
(s + 1/a! ) (s + l/a 2 ) (s + l/a n ) 


( 2 ) 


where a x , a 2 , a n , b 0 , b x , b 2 , ••*, bn_ 2 , bn-! are constants to be 

determined. Taking the partial fraction expansion of Eq. (2) gives 


H(s) = 


s + 


i/ a i 


s + 


i7 ^ + 


'n 


s + 


1/a-n 


(3) 


Then, let the input be an impulse function 

<p(t) = 7$(t ) (4) 

where s(t) is the unit impulse function and y is the magnitude of the 
impulse. The corresponding Laplace transform function of q>(t) is 

®(s) = 7- (5) 

Substituting Eqs. (3) and (5) into Eq. (l), then 

V(s) = ®(s)H(s) 

" 7 (s + 1/a! + s + l/a 2 + •** + s + l/a n )‘ ^ 

The impulse response function is 

v(t) = 7 (A 1 e" t / ai + A 2 e" t//a2 + • • • + -t/a n ) (7) 
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The impulse stimulus gives rise to an after -discharge of nystagmus 
called post-rotatory nystagmus which is the impulse response. The 
actual electrical recording of the ocular nystagmus is shown in Fig. 2 
(a) which indicates the position of the eye. Figure 2 (b) gives the 
differentiation of the trace. The upper part of (a) shows the slow- 
phase velocity which follows an exponential decay with respect to time. 
The slow-phase velocity does not decay to zero, it ceases at certain 
threshold level Vm-in . The time between the onset of input stimulus and 
the cessation of nystagmus is called post-rotatory nystagmus. 


v min 


7 (A ie “ T p/ ai + A 2 e“ T p/ a2 + 


+ A n e' T p/ an ) 


1 = A i e " T p/ a i + A a e -T p / a 2 + _ . > + A n e -Xp/a n 
y v min v min v min 


(8) 


where 


Tp = duration of post-rotatory nystagmus, and 

Vm-in = the minimum slow-phase eye velocity at the 
end of nystagmus. 

If the values of y are varied and then the values of 1 fy are 
plotted in logarithmic scale vs. each corresponding T p , the time con- 
stant of each exponent can be extracted from the curve, provided they 
are adequately separated. Practically, the input stimulations are 
pulses. It is shown (see Appendix) that the same time constants occur 
whether or not the input is an impulse or a pulse of finite width. 


b. Dynamic Characteristics of the Semicircular Canals 

(i) Transfer function of the canal - -Although overall 
transfer function discovery is very important in system analysis, 
identifying each single physical component of the system would give 
more details about the system dynamics. Fortunately, based on its 
physical structures and characteristics, the mechanics of the semi- 
circular canal can be described by a simple second-order differential 
equation . 

As early as 1931 , Steinhausen pointed out that the moment of 
inertia of the endolymph, the stiffness of the cupula and the viscous 
friction of the fluid in the canal provide this system with the proper- 
ties of a heavily damped torsion pendulum. In the idealized perfectly 
circular form of a canal (Fig. 3), the fluid ring can be considered to 
rotate around the center of the canal. Therefore, Newton's second law 
states the system dynamics as 
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Jl(<p - t) + A(cp - t) = 


(9) 


where 


0 = moment of inertia of the fluid in the canal, 

H = viscous friction at unit angular velocity, and 

A = stiffness of the cupula at unit angular displacement 

Although Eq. ( 9 ) provides the correct absolute space relations, we are 
interested in the relative angular displacement g which causes the 
neural output. Thus, we define k as cupula deflection per unit relative 
angular displacement between head and fluid, and 


g = k(cp - t) 


( 10 ) 


Substituting Eq. (10) into Eq. (9) 


0g + IIg + Ag = k©cp 


( 11 ) 


This is the linear second-order differential equation for the endolymph- 
cupula system. The corresponding transfer function of Eq. (ll) is 


S(s) _ k 

$ (s ) P"+ (n/ 0 )s + (a/©) 


( 12 ) 


A second-order equation can be considered as made up of two cascaded 
first-order lags, i.e.. 


S(s) _ k 

(s + 1/T]_ ) (s + 1/T 2 ) 


Comparison of Eqs. (12 ) and (13) yields 


T x T 2 = ©/A T!+T2=n/A 


In the case of a highly overdamped system, it turns out Ti » T 2 . 
Solving Eq. (l4) 


(13) 


( 14 ) 
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Ti ~ h/a T 2 ~ e/n 

Substituting Eq. (15 ) into Eq. (13) gives 

5(b) k 

$U7 ~ (s + A/n)(s + n/e) 

(ii) Impulse response — The constants in Eq. (16) must be 
determined in order to predict any conduct of the cupula found in 
practice. One way to solve this problem is to use a series of impulse 
stimuli with different magnitudes. The constants can be found by 
analyzing the output response. Substituting Eq. (5) into Eq. (16) 


(15) 


(16) 


H ^ ** (s + A/n)(s + n/e) 

By taking partial fraction of Eq. (17) and using H/© » a/^ gives 


(17) 


a(s) « k? ^ ( s + “ s + n/e) 


The corresponding response function is 

/ -$t 

6(t) = k 7 | (e 11 

Eq. (19) can be approximated further 

6(t) | k 7 | (e " ) (20) 



which points to a strictly exponential decay of | with time. It is 
believed that cupula displacement regulates the rate of discharge of 
the sensory cells and ampullary nerves. In turn, the neural firing 
rate is known to regulate the magnitude of responses such as nystagmus. 
Therefore, it is reasonable to use theoretical cupula displacement 
during semicircular canal stimulation as an approximate predictor of 
nystagmus. Then, by taking logarithm on both sides of Eq. (20), the 
nystagmus expires after T p sec, 
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( 21 ) 



= J log e7 + 5 l0Se 


k0 


( 22 ) 


where | m in denotes the smallest deviation of the cupula giving rise to 
a just noticeable nystagmus and Tp is called the duration of post- 
rotatory nystagmus . 

The duration Tp can be plotted in a diagram called "cupulogram" 
which was originally developed by Van Egmond, Groen, and Jongkees. The 
steepness of the nystagmus cupulogram should give the value of H/a. 

(iii) Pulse response — Theoretically, an impulse of 
acceleration requires instantaneous stopping of the skull. The term 
"impulse" here means that the braking time should be inf initesmal . This 
can not be achieved practically even using an electromagnetic brake. 

But usually a short pulse with duration less than one-tenth of the re- 
sponse time is considered as an "impulse". When the pulse width can 
not be further reduced to meet the one-tenth requirement, due to the 
limitations of the instrument, the "impulse" will become impractical. 

One way to avoid the difficulty in this time constant discovery problem 
is to let the width of all pulses be the same. 

An angular acceleration pulse with duration i and magnitude a is 
shown in Fig. 4 (a). 


cp(t) = afu(t) - u(t - t)] 


(23) 


where 


u(t - t) = 1 when 

= 0 when 

By taking Laplace transform of Eq. (23) 


(t - t) > 0 
(t - t) < 0 




( 24 ) 


Substituting Eq. (24) into (l6) 


28 




3(s) = ka(l - r^i-T — t-y 

[s(s + a/ii)(s + n/e)_ 


( 25 ) 


Taking the partial fraction expansion of Eq. (25) and approximating it 
by n/e » a/H 


H(s) 



1 - e TS \ 
s + a/II/ 


( 26 ) 


The corresponding pulse response function Fig. 4 (b) is 


|(t) = ka 







t) 


(27) 


Let t = Tp, £(T p ) = ^min, and T p > t. Equation (27) can be simplified 
as 


£min = ka | e 11 ? 


C^-x) 


( 28 ) 


The relation between T p and a is 


Tp = J log e a + fl log^g-) + “ l°SeC n T - l) (29) 

/ — T 

This equation, except for H log P \e n - 1 

A 

side, is equivalent to Eq. (22). If t , the width of the pulse, is kept 
constant, H/a should be given by plotting the cupulogram. 

Even the t's are not the same for all different magnitude pulses, 
Eq. (29) can still be modified to given the time constant, II/a, thro u g h 
an iteration process. Let 


appearing on the right hand 


C = 




T 



(30) 
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Equation (29) becomes 

(T p - t) - C = J log e a + It log e (JZ®—\ (31) 

A A X Ai nri n / 


Thus, a predictor-corrector method which is often used in numerical 
analysis to solve a differential equation can be similarly applied here 
to find the value of H/a. Let 


predictor t 

T p - t = J3 log e a + | iog e (-M_) 
A A V Almin ' 

(31a) 

corrector : 

C = | log e (e^ - 1 ) - t 

(30) 

(Tp - t) 

- C = H log e a + 5 log e f ■ 

A e A VASmin/ 

(3D 


Equation (31a) provides the starting value of H/a by graphical method 
(cupulogram) ; we then evaluate a set of C's from Eq. (30 ) and substitute 
' in Eq. (31) to generate a new value of H/a graphically. Next we evaluate 

a new set of C's based on the new H/a and again use Eq. (31) to obtain 
a next L/a. This iteration process is terminated when two successive 
iterates agree to the desired accuracy. 

(iv) Step response — The vestibular model can also be 
tested by using steps of constant acceleration and observing the latency 
between the onset of constant angular acceleration and the beginning of 
the nystagmus. Fig. 2 (b). The latency may be attributed to the fact 
that a threshold nonlinearity exists in the eye dynamics or in the C .N .S . 

The step response can be obtained by letting t approach infinity 
in Eq. (27), therefore 


S(t) = ka | 


(i 



u(t) 


(32) 


since 


lim u(t - t) = 0 

T — * oo 
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Assume |th is the corresponding threshold deviation of the cupula 
and Tp is the latency of nystagmus. 


Sth = ka | 

A 




(33) 


If a reaches ct^ which is the constant angular acceleration re- 
quired to produce a threshold deviation, the latency time, T L , will 
approach infinity. Thus, 


Hh = ka th | 


Substitute (-^h in Eq. (33) by 

Qfch n 

= 1 - e 

a 


The threshold acceleration 0%^ can be given by having Eq. (35) best 
fit the experimental data, provided JI/a is known. 

(v) The oretical value of 0, H and A— When the anatomical 
data concerning the semicircular canal are known, the constants of the 
fundamental differential equation, Eq. (ll), of the cupula- endolymph 
system may be calculated on the basis of fluid dynamics 

The moment of inertia 0 is evaluated on the assumption that the 
endolymph ring keeps the same cross section even in the utricle. That 
is, the canal extends itself through the utricle, contributing to the 
moment of inertia but not to the friction. Then, 0 (from Van Egmond) 
is expressed as 


(34) 


(35) 


0 = 2jt 2 or 2 R 3 (36) 

where 

a = density of the endolymph, 
r = cross-section radius of the canal, and 
R = radius of curvature of the canal. 
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The viscous friction n is given by 


rt = 8rtT)R 2 L 


(37) 


■where 


t) = viscosity of endolymph, and 
L = length of canal. 

The mean canal length of the rat, 1.28 itR, is approximately two-thirds 
of its circumference. Ry substituting L = 1.28 jtR into Eq. (37) 


n = 10.24 it 2 T)R 3 


(38) 


The stiffness of the cupula A is given by 


A = «pr 2 R 


(39) 


where p is the coefficient of elasticity. 

Therefore, dividing H and 0 by A and II, respectively, gives the 
value of two time constants Eq. (15) 


II 10.24 rttiR 2 
~ A “ y.r 2 


(40) 


m _ 0 _ or 2 
2 n 5.12 T) 


(41) 


RESULTS 


TIME CONSTANTS OF THE SYSTEM 

The duration of post-rotatory nystagmus for each pulse input is 
determined. Tables 1 and 2 list the means and standard deviations of 
duration for lg and 2g rats, respectively. 
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Table 1 The means and standard deviations of duration 
of post-rotatory nystagmus for lg rats 


DECELERATION 

RATE 

DURATION OF POST -ROTATORY NYSTAGMUS (sec) 

(°/sec2) 

MEAN 

STANDARD DEVIATION 

20 

4.03 

0.73 

25 

4.06 

0.72 

30 

3.41 

0.53 

40 

4.34 

1.11 

50 

4.84 

1.02 

60 

4.46 

0.97 

80 

6.16 

1.25 

100 

6.36 

2.09 

125 

7.80 

3.51 

150 

9.04 

2.40 

200 

11.50 

3.26 

250 

14.32 

3.37 

300 

15.12 

2.82 


Table 2 The means and standard deviations of duration 
of post-rotatory nystagmus for 2g rats 


DECELERATION 

RATE 

DURATION OF POST-ROTATORY NYSTAGMUS (sec) 

(°/sec 2 ) 

MEAN 

STANDARD DEVIATION 

50 

3.77 

0.57 

60 

3.85 

1.09 

80 

3.92 

0.68 

100 

4.31 

0.65 

125 

5.00 

1.07 

150 

5.04 

1.14 

200 

6.68 

1.36 

250 

8.71 

2.06 

300 

10.09 

1.40 
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The durations as a function of the inverse of deceleration rate 
are plotted in Fig. 5 (lg) and. Fig. 6 (2g), which are essentially the 
same as cupulograms. The solid curves are chosen to best-fit the data. 
If the vestibular system could be simply described by a second-order 
equation (Eq. 11), the cupulogram corresponding to it should reveal a 
straight line response (Eq. 22). However, since the response in Fig. 5 
is considered as a second-order curve by regression analysis, it sug- 
gests that the overall response includes more than one exponential term 
(Eq. 8 and A-9) . 

A logarithmic ordinate is chosen so that any single exponential 
decay should be revealed as a straight line. Although the curve in 
Fig. 5 is curved initially, it does eventually tend to a straight line, 
suggesting that this part of the response can be approximated by the 
slowest exponential term, because all faster exponentials have decayed 
to zero by this time. The slope of this straight portion gives the 
corresponding time constant . The straight part is then subtracted from 
the over-all response curve. The new curve is superimposed on Fig. 5- 
Again, as Tp increases, a single exponential is present with a shorter 
time constant. After a further repetition, a straight line is obtained, 
so that the process ends. The same technique is then applied in Fig. 6. 
Table 3 summarizes the time constants from the response curves of lg 
and 2g rats. 

A different set of deceleration rates with varying pulse widths 
were used to test lg rats. The corresponding results are summarized 
in Table 4. A method called the predictor-corrector method, which is 
described by Eqs. (30), (31), and (31a), can be applied to evaluate the 
H value of H/a. Figure 7 shows the procedures of the method and indi- 
cates the resulting II/a (9-1 sec., approximately). 

Using Eq. (4l), ®/H = ar 2 /5.12 q, with a = 1.0087 g/cm 3 , q = 0.0115 
g/cm and 7 = 0.0091 cm, gives us another time constant, T 2 = 1.42 x 10~ 3 
sec, of the semicircular canal. Since this value is too small, none of 
the suggested methods which were described previously can reveal it. 

The time constants of the semicircular canal can be summarized as 
Ti = JI/a =9*1 sec and T 2 = ®/ n = 1.42 x 10 -3 sec. 


INPUT THRESHOLD LEVEL OF THE SYSTEM 

The mean latency of nystagmus corresponding to each input deceler- 
ation is listed in Tables 5 and 6 for lg and 2g rats, respectively. 

The latencies, as a function of deceleration, are shown in Figs. 8 
and 9. By using Eq. (35) with known K/a = 9*05 for lg rats and H/a = 8.42 
for the 2g rats, the value can be adjusted until the equation fits the 
experimental data. The solid curve in Fig. 8 gives at = 5°/sec 2 for 
1G rats and the one in Fig. 9 shows a higher input threshold level 
(11° /sec 2 ) for the 2g rats. 
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INVERSE OF DECELERATION RATE IN SEC 2 /. 



Table 3 The time constants of lg and 2g rats 



TIME 

CONSTANT (sec) 



a l 

a 2 

a 3 

1G 

9.05 

1.30 

0.39 

2G 

8.41 


0.30 


Table 4 Mean duration of post-rotatory 
nystagmus (Tp) for lg rats 


DECELERATION 

RATE 

DECELERATION 

TIME 

MEM DURATION OF 
POST -ROTATORY NYSTAGMUS 

(°/scc 2 ) 

(sec) 

T 

(sec) 

T P 

T P ' T 

30 

20.0 

21.43 

1.43 

40 

15.0 

17.78 

2.78 

50 

12.0 

15.91 

3.91 

60 

10.0 

15.50 

5.50 

80 

7.5 

14.09 

6.59 

100 

6.0 

13.71 

7.71 

125 

4.8 

13.09 

8.29 

150 

4.0 

13.05 

9.05 

130 

3.33 

13.01 

9.68 

210 

2.86 

12.67 

9.81 
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Table 5 The mean and standard deviation of latency 
of nystagmus for lg rats 


DECELERATION RATE 

LATENCY OF NYSTAGMUS (sec) 

(°/sec^ ) 

MEAN 

STANDARD DEVIATION 

20 

2.96 

2.16 

25 

1.20 

0.35 

30 

1.39 

0.64 

40 

1.00 

0.43 

50 

0.75 

0.14 

60 

0.77 

0.33 

80 

0.63 

0.25 

100 

0.56 

0.21 

125 

0.40 

0.09 

150 

0.37 

0.11 

200 

0.27 

0.06 

250 

0.25 

0.06 

300 

0.21 

0.08 
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Table 6 The mean and standard deviation of latency 
of nystagmus for 2g rats 


DECELERATION RATE 

LATENCY OF 

NYSTAGMUS (sec) 

(°/sec2) 

MEAN 

STANDARD DEVIATION 

20 

12.82 

5.93 

30 

2.90 

1.36 

40 

2.72 

1.58 

50 

1.35 

0.17 

60 

1.64 

0.62 

80 

1.20 

0.50 

100 

0.98 

0.38 

125 

0.95 

0.36 

150 

0.80 

0.31 

200 

0.56 

0.16 

250 

0.41 

0.17 

300 

0.39 

0.15 
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Fig. 8 The latency of nystagmus following the onset 
of constant angular deceleration for lg rats . 
The threshold deceleration is 5°/sec 2 . 
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LATENCY IN SECONDS 



60 120 180 240 300 

RATE OF DECELERATION IN °/SEC 2 


Fig. *9 The latency of nystagmus for 2g rats. The 
threshold deceleration indicates a higher 
threshold level of 11° /sec 2 . 
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DISCUSSION AND CONCLUSION 


One of the primary purposes of this study is to evaluate pulse 
wave forms of angular decelerations and the response they elicit as a 
means of discovering the time constants of system in vestibular nystagmus 
response. Pulse waveforms were selected as a test stimuli because they 
eliminated the difficulties due to the limitations of the stimulator 
and provided an adequate method to test the rats . 

However, the main goal is to find the response difference between 
rats chronically exposed to a lg and 2g environment on a basis of 
biocybemetic models of the horizontal vestibulo- ocular reflex arc. 

Before presenting the models, the time constants in Table 3 should 
first be related to the physical model. 

Deceleration rates with varying pulse width were selected as 
stimuli. The advantage of using longer widths for low deceleration 
rates (Table 4) is to minimize the effect of exponential terms with 
shorter time constants on the overall response. Therefore, the system 
tends to respond as a first-order system under these special stimula- 
tions. Although the cupulogram (Fig. 7) showed a straight line response, 
corrections are still needed to evaluate the actual value of n/A. 

The value of II/a, 9-1 sec, is close to both a^s. We can assume a x 
is attributed by the semicircular canal. H/a is determined by the 
canal's physical dimensions Eq. (40) which are not changed during 
chronic centrifugation. Assume the physical properties of canal, p, q 
and cr, are not changed. The value of h/a in either state should be 
equal to each other. Again the values of a x shows the validity of this 
assumption. 

Comparing the time constants in Table 4 with the block diagrams in 
Fig. 6 (a) leads to two assumptions . First, a 3 may be related to eye 
dynamics which is not affected by the chronic 2g state. Ten Kate sug- 
gested a value between 0.25 and 3-97 sec for the time constant of pike's 
eye dynamics and a value of 0.32 sec for the cat. This may further 
support the hypothesis. 

a 2 may be contributed by the utriculus via the C.N.S. Since linear 
acceleration influences the responses of the semicircular canals, a 
prolonged 2g state may interfere in some way with the vestibular re- 
sponse. The resulting effect could be the reason for the missing a 2 
in the response of 2g rats. 

The forces produced by gravitation, linear acceleration, and cen- 
trifugation are identical with respect to their effect on organisms . 
However, centrifugation differs from the other in that it involves a 
rotation. The coriolis acceleration produced by body motions in a 
rotating environment may habituate the nystagmus response and lead to 
the absence of a 2 . This will need to be examined further. 
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The time constants of the semicircular canal were found as Tx = 9*1 
sec and T 2 = 1.42 x 10“ 3 sec. The inverse of T x and T 2 are corresponding 
to the lower and higher 3-hB frequencies of the canal in frequency domain. 
Therefore, the bandwidth, i.e., the range of frequency between higher and 
lower 3-dB frequency, is roughly from 0.11 rad/sec to 705 rad/sec. Mayne 
stated the bandwidth of the canal appears to cover the range of body 
movement. Due to the small size and light weight of the rat, its band- 
width should include a range higher than that of a human being (from 
0.04 rad/sec to 300 rad/sec). 

Summarizing the discussion and the results in the previous section, 
the mathematical models of lg and 2g for acceleration input are shown 
in Fig. 10. The value of 9 is chosen as a nominal value for the semi- 
circular canal. K 1} K 2 and K 3 are gains for each block. The model for 
canal dynamics follows Eq. (12). The one suggested for eye dynamics 
describes the damping of the eye muscles, however. To is only an assump- 
tion. Since the denominator is first order, a first-order term may or 
may not exist in the numerator. 

The threshold which only represents an over- all response of the 
system is not consistent with any physical structure in the horizontal 
vestibulo-ocular reflex arc. 

Numerous investigators have suggested that the responses of the 
semicircular canals may be described by a heavily over damped torsion 
pendulum. This means a nystagmus cupulogram should reveal a straight 
line response rather than a curve response which was shown in Figs. 6 
and 7. It appears that this assumption is too simple and that the 
deviations from the presupposed mechanical laws point out either other 
physiological structures may play an important role in vestibular 
nystagmus response or the semicircular canal itself can not be approxi- 
mated by a simple linear second-order model. The semicircular canal 
structure has its nonlinearities in the sense of mechanics. The spring 
constant of the cupula may not be constant even for small deviations 
from the resting position; it may be a function of its own position. 

The damping factor is not linear in its nature and may be in quadratic 
form or have directional preference. 

The difference between the responses of lg and 2g rats may also be 
explained by the nonlinear characteristics of the semicircular canals. 

It is possible that chronic centrifugation may have an effect on the 
vestibular organs . Anatomical studies on rats exposed to chronic 
acceleration have shewn considerable changes in body size and in tissue 
composition. These changes may result from the increased apparent 
density of body tissues and fluids induced by the greater accelerative 
force. Similarly induced changes in chemical composition of the laby- 
rinth may have resulted from chronic hypergravity. Therefore the 
density and viscosity of endolymph or the spring coefficient of the 
cupula may differ from those of lg rats. These differences, if they 
exist, could alter the functional characteristics of the semicircular 
canal. 
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Fig. 10 (a) The mathematical model of lg for acceleration input, (b) The 

model of 2g with a missing utricular term and a higher threshold. 












APPENDIX 


An angular acceleration pulse of duration r and magnitude a can be 
expressed as 


<p(t ) = a[u(t) - u(t - t)] 


(A-l) 


The corresponding Laplace transform is 


$(s) = TS - 


(A-2) 


Substituting Eq. (A-2) in 


V(s) = $(s) 


Vi ® 11 " 1 + 


+ b^s + b c 


_( s + 1/ai )(s + l/a 2 ) ••• (s + l/a n )_ 


(A-3) 


gives 


tr , ^ /, — ts \ l~ b n s 11-1 + ••• + b s + b 

( s ) — Q '( e ^_s(s + l/ai)(s + l/aa) • (s + l/a n )_ 


(A-4) 


Take partial fraction of Eq. (A-4) 


V(s) = cc(l - e TS )^— + — + ... + %— ) 

\ S S + 1/9*1 s + 1/ 9*2 s + l/^n/ 

(A-5) 

The corresponding impulse response function is 

v(t) = ako [u(t) - u (t - t)] + aAl £e u(t) - e ( t -' r V ai u (t _ T )J 
+ aAa [e^^ 8,2 u(t) - e “(’ fc-T )/ a 2 _ T jJ + ... 

+ akn [e-^ u(t) - e -(t-T) / a n u(t . T) ] (A-6) 
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Let Tp be the duration of post-rotatory nystagmus and Vm-in be the 
minimum slow-phase eye velocity. Then substitute Tp and v-min in 
Eq. (A- 6 ) 


v min = aAo [u(T p ) - u(T p - t)] + aAi [e Tp//ai u(T p ) - e ^ Tp u (T p _ x 

+ ak 2 [e“ Tp ^ ai u(T p ) - e ”^ T P"' r V a 2 u ( Tp _ 

+ ... AA[e' Tp / an u(T p ) - e ~ (T P~ T V a n u(T p - t)] (A-7) 

Since Tp > t, Eq. (A-7) becomes 

Vmjn = o:[ai(i - e T / ai ) e _T p/ ai + Aa(l - e T / a 2 ) e ~ T p / a2 + • •• 

+ AA(l - e T / %1 )e' T p/ a ^] (A-8) 

Divide both sides of Eq. (A-8) by 


- = - (l 

a Vmin 


e Va 1 ) e -Tp/a 1 + _A|_ ^ _ e Va 2 ^-Tp/a 2 + _ 
v min 


+ -is. (l - e T /%>) e -’'p/an 

Vmin v 


(A-9) 


The relation between a and Tp is equivalent to that of 7 and T p in 
Eq. ( 8 ) except the constant terms (l - e T / ai ), (l - e T / a2 ), •••, 

(l - e ^ &n ) existed in Eq. (A-9). Therefore, the same method can be 
applied to find the time constants in pulse response. 
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